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Preface 


A means  of  passively  tracking  a target  in  range  has  been  sought 


since  the  beginning  of  tracking  systems.  Recent  develojjments  in  solid 
state  electro-optic  detectors  have  made  very  accurate  angle  tracking 


possible.  At  the  same  time,  solid  state  technology  has  also  made  it 
possible  to  process  video  data  very  rapidly  using  small  devices.  The 
result  is  a small,  light-weight,  high  resolution  sensor  capable  of 
following  a target  at  high  rates.  The  impact  of  aircraft  vibrations 
on  those  sensors  is  not  well  known.  The  purpose  of  this  investigation 
is  to  determine  the  effects  that  the  vibration  environment  will  have 
on  the  functioning  of  these  devices. 
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Abstract 

The  stereometric  ranqe  t indinq  technique  is;  used  as  the  basis 
oi  an  airborne  ivissivo  tracker.  An  analysis  of  the  method  reveals 
that  anqul.u  accuracy  is  the  most  critical  element  in  this  technique. 
Uanae  accuracy  can  be  improved  by  increasinq  baseline  separation. 
However,  t tie  vibration  environment  of  the  aircraft  becomes  worse  as 
tlu'  baseline  is  increased.  The  results  ot  this  investigation  indicate 
that  these  vibrations  can  counter  most  of  the  qain  in  tesolution.  The 
invest  iq.it  ion  is  hindered  by  a lack  of  vibrational  data,  especially 
anqulat  vibration  data.  However,  the  invest  in.it  ion  establishes  what 
data  is  needed,  and  how  such  data  can  be  applied,  once  obtained,  to 
estimate  I'ortormance  potential  of  the  stereometric  ranqinq  method. 


i .troduction 
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Backg  rou  mi 

F' ight  or  aircraft  perform  a number  of  different  missions.  These 
missions  can  be  divided  into  air-to-air  and  air-to-ground  missions. 

In  air-to-air  missions,  the  target  is  an  airborne  object  such  as  another 
aircraft  or  missile.  Intercepts,  escort,  and  close-in  aerial  combat 
are  examples  of  air-to-air  missions.  In  air-to-ground  missions,  the 
fighter’s  target  is  on  the  ground.  Such  targets  include  enemy  tanks, 
bridges,  air  fields  and  logistics  areas.  Air-to-ground  missions  include 
close-air  support,  interdiction  and  airfield  attack. 

Whether  the  target  is  on  the  ground  or  in  the  air,  the  aircrew's 
mission  is  to  find  it  and  put  ordnance  on  it  when  required.  In  order 
to  successfully  complete  the  mission  without  unnecessary  risks  to 
aircraft  and  crews,  sophisticated  tracking,  fire  control  and  weapons 
release  systems  are  installed  on  USAF  fighters. 

Tracking  systems  are  used  to  follow  the  motion  of  the  target 
relative  to  the  fighter.  Information  from  the  tracking  systems  is 
processed  to  determine  target  range,  range  rate,  line-of-sight  (LOS) 
angle  and  angular  rate.  Fire  control  systems  present  this  information 
to  the  crow  in  a suitable  format.  Weapon  release  systems  use  tracking 
information  along  with  fighter  state  data  to  compute  weapon  release 
or  to  launch  signals. 

Tracking  systems  can  be  classified  as  passive  or  active.  Active 
tracking  systems  are  those  winch  direct  energy  at  the  target  and  track 
reflected  energy.  Radars  and  laser  designators  are  examples  of  such 
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systems.  Passive  tracking  systems  are  those  which  track  enerqy 
emitted  by  the  target.  Heat  seeking  devices  and  electro-optic 
trackers  are  examples  of  passive  trackers. 

At  present,  fighters  must  use  active  trackers  to  measure  target 
range  and  range  rate  in  real  time.  However,  active  trackers  alert 
the  target  to  the  presence  of  the  fighter.  This  gives  the  enemy  forces 
time  to  take  evasive  maneuvers,  use  electronic  countermeasures,  fire 
homing  missiles  back  at  the  fighter  or  take  cover.  The  element  of 
surprise  is  lost.  Not  only  does  this  decrease  the  chances  of  destroying 
the  target,  it  also  increases  the  hazards  to  the  attacking  force. 

For  example,  according  to  Reference  1,  the  McDonnell  Douglas 
F-4  rhantom  has  only  limited  capability  to  attack  a MIG  25  Foxbat  at 
80,000  feet  and  Mach  2.0+.  This  would  surely  be  negated  if  the  F-4 
were  to  use  its  radar,  an  active  system,  and  alert  the  Foxbat  that  it 
was  under  attack.  But  the  F-4  must  use  its  radar  to  detect  and  track 
the  Foxbat.  For  this  example,  a passive  tracking  system  which  has  a 
range  tracking  capability  would  increase  the  success  of  the  mission. 

It  would  be  possible  to  use  the  tracking  information  to  intercept  the 
target  at  optimum  missile  launch  parameters  before  turning  on  the  radar 
for  missile  guidance.  The  first  indication  of  the  attack  would  be 
the  missile  launch. 

Another  example  of  the  advantage  of  a passive  tracker  is  the  case  of 
a single  fighter  attacking  a ground  target  deep  in  enemy  territory.  For 
mission  success  and  fighter  survival,  such  missions  rely  on  one  surprise 
pass  at  the  target.  This  means  that  the  fighter  should  not  emit  any 
energy  as  it  navigates  to  the  target.  Accuracy  of  the  bomb  run  is 
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important  since  the  aircrew  can  afford  to  make  only  one  pass.  On  the 


bomb  run,  the  crew  would  like  to  go  in  as  fast  as  possible,  fly  an 
unpredictable  flight  path  and  keep  look-out  for  enemy  surface-to-air 
missiles  (SAM) . This  means  that  the  actual  bomb  release  would  best 
be  done  by  an  automatic  weapons  release  system.  But  such  systems 
often  rely  on  target  range  information  from  an  active  tracking  system. 
Such  a system  is  vulnerable  to  countermeasures.  A passive  tracking 
system  is  less  vulnerable  and  almost  impossible  to  detect. 

Present  passive  tracking  systems  in  fighters  are  LOS  trackers 
only.  What  is  needed  is  a passive  system  which  is  also  capable  of  range 
tracking.  An  extension  of  the  LOS  tracker  would  be  a stereometric 
system  composed  of  two  such  trackers  separated  by  a baseline.  It  will 
be  shown  later  that  this  system  requires  very  accurate  LOS  performance 
from  the  trackers. 

Purpose 

It  is  the  purpose  of  this  report  to  perform  a parametric  trade-off 
analysis  of  a passive  stereometric  range  finding  system.  Sensor  angular 
resolution,  aircraft  vibration  and  flexures,  and  sensor  baseline  will 
be  analyzed  to  determine  their  effect  on  the  range  finding  accuracy  of 
this  system. 

This  analysis  will  be  applied  to  a specific  fighter  aircraft  to 
determine  sensor  specifications  required  to  achieve  practical  range 
accuracies  for  various  air-to-air  and  air-to-ground  missions. 

Approach 

This  analysis  will  start  with  a discussion  of  the  theory  of 
stereometric  ranging.  An  error  analysis  will  be  performed  to  disclose 
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the  most  important  par.uneters  for  range  accuracy. 

Some  criteria  foi  tracker  performance  are  assessed,  starting 
with  contract  specif icat ions  for  research  on  an  advanced  stereometric 
tracker.  The  criteria  are  expanded  to  meet  the  air-to-air  and  air-to- 
ground  mission  requirements. 

Angular  resolution  of  a focal  plane  sensor  is  investigated  as  well 
as  the  means  used  to  track  the  target. 

t 

In  order  to  increase  maximum  range  capability  of  the  tracker, 
the  effect  of  placing  two  trackers  at  remote  locations  on  the  aircraft 
is  studied.  Some  typical  aircraft  dimensions  are  given. 

The  effect  of  aircraft  vibration  and  flexures  are  presented. 
Analytical  equations  are  formed  whenever  possible.  The  expressions 
derived  from  this  analysis  are  used  to  evaluate  sensor  placement  on 
a typical  fighter. 

The  last  section  contains  the  conclusions  from  the  analysis  and 
appropriate  recommendations  are  made  for  the  development  of  the 
stereometric  tracker. 


•1 


II.  St  el  outlet  I if  Kanqe  Method 


lilt  I odltf  l loll 

St  oi oomet  i if  ramie  theory  ia  moat  easily  presented  in  qeometr  if 
lotma.  one  hum  mu  i ea  tlio  lino  ot  aiqht  (U'S)  amilo  l o ,m  obioct  trout 
two  posit  Iona  a lomj  .1  baseline,  In  l'iqure  1,  t ho  l .OS  am)  l o is  measured 
with  Uio  hol|>  ot  a roteieneo  lino,  Kl.  . Thin  lino  ia  perpendicular 
to  t ho  Imho  amt  in  I ho  piano  formed  t>y  the  baseline  and  tarqet, 

I’aiallaK  ia  measured  hy  aliiinimi  t ho  tip  ot  the  reference  lino  and  t ho 


t arqot  . 


‘I' ho  difference  in  the  1.0S  ami  lea  ia  a . In  the  special  case  when 
the  tarqet  l .OS  ia  pai  allot  to  the  reference  lino  (t'iq.  1),  it  ia  oaay 
to  see  that 


t all 


.i 


II 

II 


l’ 

III. 


(1) 


whoi  o 


.i  - ann to  formed  hy  both  l.OS's  at  the  tarqet 

ti  - baseline  separation 

K.,  - ranuo  to  the  tarqet  from  position  •! 

I'  - parallax  in  shjhtinq  amito  duo  to  a shift  in  viewinq  position 

Kl.  - lonqth  ot  reference  lino 


•> 
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A 1 so , with  the  hel|>  of  Figure  1 


is 


1 1 


a in  a - tan  u 


t lion 


is 


R 


1 


B = V_ 
R Rl, 


(2) 


where 


R 


1 


range  to  target  from  position  1. 


The  approximation  is  valid  as  long  as 


sin  a 

tan  a = 3 sin  a 

cos  a 


or 


cos  u 
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(3) 


or 
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will' I e 


R i s the 


For  the  purpose  of  this  analysis,  R^  “ R,  ■ R , 
ranqe  to  the  tarqot  item  either  sensor.  Equation  2 will  bo  written  as 


B I' 

— m — 

R F 


>1 


(4) 


Error  Analysis 

This  paper  will  analyze  the  ranqe  trackinq  capability  of  a 
stereometi ic  system  mounted  on  a I iqhtet  aircraft.  Such  a system  will 
be  composed  of  two  sensors  capable  of  measurin'!  imaqe  shift  due  to 
baseline  separation  (see  Fiqure  2).  It  is  important  to  know  the 
par.imeters  in  Equation  (•!)  which  will  limit  maximum  ranqe  and  ranqe 
accuracy  of  the  system, 
bet 


AR 

a 

ranqe  accuracy  desired 

AH 

aecu i acy 

of 

base) ine  measurement 

AF 

rs 

accuracy 

of 

focal  lcnqth  measurement 

AH 

•a 

accuracy 

of 

pai al lax  measurement 

Substitutin'!  these  values  into  Equation  (4)  qives 


K i A K 
B t AH 


V f AF 

p r ae 


K * AR 


EH  t A EH  1 _E AH  * ARAB 

(r  + ah) 


H 


POINT  SOliKCK 
^ > OF  ENERGY 

I i 


n 

■ ’ 

/T*  °» 

vb 


THE 


FUUIUI 


Kami*'  Mi>.i:iurt'im>nt  by  Stereomoti  ic  Method 


Assuming  At’  <<  r , 


l 


1 


P (l  ♦ A-;  > P 


1 

cr  — - 

p 


AP  /AP\-'  /ap\ 

- T * (Tj  - \T) 


1 - 


AP 


Substitution  into  the  previous  equation  qives 


K t AR  ■ | KB  + AFB  + FAB  r AKAB] 


s (■  - f) 


FB  AFB  FAB  AFAB 

P P P P 


- AP 


PB  AP 

' p 

P^ 


AFB  FAB  AFAB 
P * P + P 


Assuminq  AF  <'  F and  AB  <<  b , the  last  term  in  this  equation 


AP 

P 


AFB  FAB  AFAB 

P + P P 


is  insiqnif ieant  compared  to  the  other  terms  in  the  equation. 
Thus, 


FB  . AFB  FAB  (-)FBAP 

* ♦ AR  - -p-  ♦ -p“  + -pT  + — 


p 
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Since 


FB 


AR 


AFB  FAB  FBAP 

P P + „2 


AR 

R 


AF  AB  AP 

F B P 


In  words  then,  the  sum  of  the  parameters  expressed  as  ratios  is  equal 
to  the  range  error  expressed  as  a ratio. 

This  analysis  is  concerned  with  angular  resolution,  a , and 

baseline  separation  on  ranging  accuracy.  Therefore,  the  effects  of 

AB  , AF  . , , , . . . 

— and  — on  range  accuracy  will  be  ignored.  As  a word  of  caution, 

B F 

however,  the  error  in  F and  B cannot  be  ignored  when  specifying 
the  tolerances  for  sensor  optics  and  baseline  separation.  It  will 
always  be  at  least 

It  appears  that  range  accuracy  is  not  affected  by  the  error  in 
parallax  any  more  than  it  is  affected  by  error  in  baseline.  However, 
if  baseline  and  parallax  are  expressed  in  terms  of  range,  then 


AF  All  < AR 
F + B ~ R 


AR  _ AB 
R B 


(5) 
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and 


AH 

R 


AT 

P 


1 

since  — 


R 

KB 


AR 


R2 


P 

FB 


(6) 


Tlie  quantity  AR  increases  linearly  with  R , once  AU  and  B are 
set.  However,  AR  increases  with  the  square  of  R if  AP , F and 
B are  constant. 

AR 

If  the  ratio  — is  made  a constant  for  all  ranges,  AB  is  also 

AR 

a constant.  However,  keeping  — constant  for  all  ranges  does  not 

R 

make  AP  a constant.  It  still  varies  linearly  with  range  and  it  will 
be  the  most  critical  parameter  at  long  ranges.  If  a particular  system 
is  capable  of  measuring  ixirallax  to  AP  and  the  system  is  required  to 
measure  range  to  a qiven  accuracy,  expressed  as  AR/R  , then  the  maximum 
range  of  the  system  is  set  by  Equation  (6) 


R 

max 


(6a) 


Summary 

This  section  examined  the  stereometric  range  theory.  Several 
assumptions  were  necessary  to  maintain  a simplified  form  of  the  ranging 
equation. 

An  error  analysis  was  performed  on  the  ranging  equation.  In 
order  to  extend  the  range  capability  of  system,  it  appears  that  AP 
will  prove  to  be  the  limiting  term. 


III.  Tracker  Requirements 

This  section  will  list  performance  criteria  for  the  tracker. 

These  criteria  are  based  on  a statement  of  work  for  Contract  F3361S- 
78-C-1S62  (Ref.  1).  The  criteria  stated  in  that  document  are  for  close- 
in  aerial  gunfire.  This  environment  makes  the  greatest  demand  on  the 
tracker  in  terms  of  angle  rates.  Maximum  range  of  current  air-to-air 
intercept  missiles  are  used  to  estimate  desired  maximum  range  capability 
of  the  tracker  (Refs.  2,  3,  -1).  The  desired  field  of  regard  is  a 
combination  of  the  field  of  regard  required  in  the  contract  and  the  field 
of  regard  for  the  Pave  Tack  system  quote  in  Reference  5,  page  57. 

Table  I shows  the  parameter  requirements  from  Reference  1. 

Table  I.  I.ist  of  CAI  Criteria  (Ref.  l:A6-7) 


Parameter 

Magnitude 

Accuracy 

Field  of  Regard 

120°  Cone 

4 1 milliradian 

Angle  Rate 

Maximum  of  60°/sec 

4 4 milliradians/sec 

Maximum  Range 

4000  ft 

4 + 50  ft 

Range  Rate 

+500  to  -1600  ft/sec 

4 ± 50  ft/sec 

The  field  of  regard  parameter  of  a 120°  cone  is  referenced  to  the 
aircraft  fuselage  longitudinal  axis  (see  Figure  3) . The  article  on 
Pave  Tack  states  that  that  system  has  a full  hemispherical  field  of 
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FIGURE  3.  120°  Cone 


regard  below  the  aircraft.  The  field  of  regard  requirement,  then,  will 


expand  to  include  the  lower  hemisphere  also.  The  maximum  range  require- 
ment has  little  use  outside  of  an  air-to-air  gun  tracking  mission. 

Table  II  gives  maximum  ranges  of  some  current  air-to-air  guided  missiles 
from  unclassified  sources. 

Table  II.  Missile  Ranges 


Name 

Designation 

Maximum  Range 

Falcon 

AIM-4A 

5+NM  (Ref  4:82) 

Sidewinder 

AIM-90 

10NM  (Ref  4:114) 

Sparrow 

AIM-7F 

24+NM  (Ref  3:82) 

Phoenix 

AIM- 54 

89NM  (Ref  2:56) 

To  be  compatible  with  these  missiles  and  future  missiles,  the  tracker 
should  have  a maximum  range  of  89  NM.  It  is  often  desirable  to  launch 
a missile  at  a range  which  will  give  it  the  greatest  probability  for 
a kill.  If  this  optimum  range  is  assumed  to  be  half  the  maximum  range, 
then  tracking  beyond  that  range  might  be  considered  unnecessary. 

However,  a fighter  often  wishes  to  improve  its  position  with  respect 
to  the  tarqet  prior  to  missile  launch.  To  do  so  requires  time  and 
distance.  Tf  the  fighter  and  its  target  are  closing  at  1000  knots, 
they  close  16-2/3  miles  every  minute.  Therefore,  if  optimum  range  for 
the  AIM-54  is  89NM/2  and  if  the  fighter  requires  about  two  minutes  prior 
to  launch  for  maneuvering,  maximum  range  capability  of  the  tracking 
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system  is  2(16  2/3) NM  + 44  1/2NM  = 77  NM  . The  AIM-54  is  a long 

range  missile.  Maximum  tracking  for  the  other  missiles  varies  from 
45  NM  to  35  NM.  For  the  purposes  of  this  thesis,  the  maximum  tracking 
range  of  the  system  will  be  set  at  89  NM.  This  range  could  easily  be 
twice  what  is  actually  required  if  the  fighter  is  not  compatible  with 
the  AIM-54,  or  if  missile  launch  range  is  restricted  to  even  shorter 
ranges  due  to  other  factors.  Maximum  tracking  capability  has  been 
assessed  using  air-to-air  missiles  because  air-to-ground  ordnance 
typically  lias  shorter  range.  If  the  range  accuracy  of  50  ft  at 
4000  ft,  as  stated  in  the  contract,  is  extended  to  89  NM,  range 
accuracy  would  be  1.1  NM. 

The  -1600  ft/sec  range  rate  is  not  sufficient  for  two  supersonic 
fighters  closing  head-on  (Reference  6).  In  the  reference  cited,  two 
fighters  were  closing  at  Mach  4.0.  This  converts  to  a closing  rate  of 
up  to  -4350  ft/sec,  depending  on  precise  atmospheric  conditions. 
Maintaining  the  range  rate  accuracy  ratio,  the  accuracy  requirement 
is  136  ft/sec. 

Both  the  range  and  range  rate  accuracy  requirements  in  Table  I 
were  to  predict  the  path  of  a bullet  fired  from  a fighter.  The 
predicted  bullet  trajectory  is  used  to  generate  an  aiming  reference 
for  the  pilot.  This  requirment  may  be  more  restrictive  than  required 
for  a guided  missile.  Table  III  lists  revised  parameters  for  the 
stereometric  tracker. 


Table  III.  Revised  Tracker  Requirements 


Parameter 

Magnitude 

Accuracy 

Field  of  Regard 

Angle  Rate 

120°  Cone  forward 
plus  Hemisphere 
below  aircraft 

4 1 milliradian 

Angle  Rate 

60°/sec  (Max) 

4.  4 milliradians/sec 

Maximum  Range 

89NM 

4 ±1.1  NM 

Range  Rate 

500  to  -4350  ft/sec 

41136  ft/sec 

IV. 


Syst cm  Align lai 

Inti  oduct  ion 

Maximum  range  capability  and  range  accuracy  depend  on  bow  well 
the  angle  a can  be  measured.  Angle  u is  tlie  LOS  difference  between 
the  two  sensors  (Figure  4). 


‘Not 


a 


- 0. 


(6.1) 


where 

o. 

35 

LOS  angle 

measured 

1 

normal  to 

the  base 

by  sensor  1 relative  to  a vector 


0 


2 


LOS  angle  measured  by  sensor  2 relative  to  the  same 
reference 


If  0 and  0 are  measured  to  an  accuracy  of  AO,  and  AO. 

in  ■ I 


respectively,  then 


Act  = fa*  + A022  (7) 

where  Au  = accuracy  to  which  a is  known. 

Ideally,  sensors  1 and  2 are  manufactured  such  that  AO  = A0^  = A0 
and  Act  = /l  AO 

Sen s or  Resol i it  ion 

The  angular  resolution  of  the  sensor  is  determined  by  the  detector's 
smallest  resolution  element  (fief.  7:1),  For  a sensor  composed  of  a lens 
which  focuses  an  image  on  an  array  of  photosensitive  elements,  the  limit 
of  resolution  is  the  area  of  one  photosensitive  element  divided  by  the 
focal  length  squared.  The  units  arc  steradians.  If  the  image  of  the 
target  if:  smaller  than  one  resolution  element,  ,tj,'e ^response  of  that 
element  is  the  same  as  if  the  entire  element  were  illuminated  by  the 
image.  For  a stereo  range  finder,  the  angular  resolution  in  the  direction 
parallel  to  the  base  is  of  interest.  In  this  case,  it  is  meaningful  to 
describe  angular  resolution  as  the  dimension  of  a photosensitive  element 
parallel  to  the  base  (see  Figure  5) . 
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FIGURE  5.  Angular  Resolution 


Then 


. "I  P P 

AO  = sin  — — 

r r 

where 

AO  - is  measured  in  milliradians 
p = width  of  photosensitive  element 

K = focal  length  of  optical  system 

CAI  Proposal 

Tlie  manner  in  which  u is  measured  can  be  improved  using  a 
proposal  by  CAI,  Division  of  Bourns,  Inc.  (Reference  8).  In  this 
proposal,  the  image  in  one  sensor  is  compared  directly  with  the  image 
in  tlie  other  sensor  by  means  of  a correlation  algorithm.  The  resulting 
accuracy  in  n is  a fraction  of  A6  . The  following  is  a description 
of  that  proposal. 

The  passive  tracker  consists  of  three  sensor  systems  in  one 
housing  (see  Figure  6) . The  stereo  base  is  four  inches.  The  middle 
lens  focuses  the  target  on  detector  number  2.  All  three  detectors  are 
charge  coupled  devices  (CCD).  These  detectors  are  in  the  focal  plane 
of  the  lenses.  The  detectors  record  the  target  image  over  one 
integration  period  (Ref.  9:58)  and  produce  a video  signal.  The  video 
signal  from  the  center  detector  is  stored.  The  next  signal  is  compared 
to  the  stored  signal  by  means  of  a correlation  algorithm.  The 
misregistration  in  successive  image  signals  computed  using  this 
correlation  algorithm  is  used  to  produce  correction  signals  to  servo- 
mechanism controlling  the  pointing  direction  of  the  entire  housing 
(Ref  8:2.6). 
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The  rnnae  sensors  me  cimiposed  ol  ramie  lenses  and  CVP  detector 
,u  i .tys.  The  video  sianaLs  aonornted  by  detector  1 and  i are  used  to 
compute  inuuje  rnistoa  ist  i at  ion  between  images  I i om  two  sonsoi  s at  t lie 
same  time  (Rot  H:„'.n).  Since  they  are  separated  by  a baseline,  target 
I an>ie  can  be  computed  by 


K 


I1 

Kit 


whet  e K im.iao  m t si ea i st rat  ion  between  sensor  l and  i. 

iYl>  1'etectoi  s.  The  detectors  used  tor  the  ramie  sensors  are 
Fairchild  iVP  .'ll  .lira  arrays.  They  are  metal-ox ide-semiconductor 
(Mi vs)  structures  (Rot.  -li.ltl).  The  semiconductor  is  di'ped  silicon  so 
that  tin'  defcctoi  lias  .-a  I icon  j esj'onsiv  i t y (Ret.  H:»'.l-1).  The  array 
has  .’-bl  elements  by  p>0  elements.  The  active  area  ot  each  element  is 
IS  in  i ci  omet  ers  vertically  by  11  micrometers  hoi  i r.onta  1 ly  . lively  othei 
column  is  a vettical  analoa  iciintei  u.sed  to  t i uuspot  t charae  packaaes 
(see  Kiaure  7).  As  stated  before,  the  iniaae  i'i  t ocust'd  on  this  array 
by  a lens  system.  Kneh  photosensitive  element  converts  t he  inciMiiina 
lialit  enemy  into  pilot .o  election-hole  pairs.  This  process  is  called 
detection.  All  the  chaiae  carriers,  eithet  the  electrons  or  holes,  are 
accumulated  at  each  photosensitive  element  and  are  held  there  by  a 
potential  well.  They  are  accumulated  for  a period  of  2.2‘j  ms,  called 
the  int  miration  period.  At  the  end  of  the  integration  jvriod,  all  the 
packages  of  charge  are  moved  into  the  vertical  registers  .it  one  time. 
Then  the  packaaes  are  moved  one  step  up  to  the  horizontal  tea  ist or.  An 


FIGURE  7.  Diagram  of  Fairchild  211  CCD  (Ref.  1:2.15) 


entire  horizontal  line  of  charge  packages  is  transported  through  a 
preamplifier  which  sends  the  video  signal  to  an  amplifier  off  the  CCD 
chip.  Then  the  vertical  registers  move  all  the  charge  packages  up  one 
more  increment  so  that  the  next  line  of  charge  packages  is  in  the 
horizontal  register.  The  process  continues  until  all  the  charge 
packages  have  been  read.  In  effect,  the  image  is  broken  down  into 
elements,  then  these  elements  are  read  line  by  line  in  the  horizontal 
register  (Ref.  9:58).  The  time  it  takes  to  read  all  the  charge 
packages  is  referred  to  as  the  read-out  period.  It  does  not  exceed  the 
integration  period. 

The  video  signal  from  the  CCD  chip  is  further  amplified  and  then 
sent  to  an  analog-to-digital  converter.  The  digitized  video  signal 
is  then  sent  to  the  correlation  processor.  Here,  the  signal  is  compared 
with  a previously  stored  image  signal  by  a correlation  algoritlun  and 
the  image  shift  is  found. 

Correlation  Algorithm.  The  correlation  function  for  comparing 
two  scenes  is 


1^ 

A 


P 


(x,y)  Q(x  + Ax,  y + Ay)  dxdy 


where 


P (x,y)  = a function  representing  the  target 

image  over  the  detector  surface 

Q (x  + Ax,  y + Ay)=  a function  representing  the  target 

image  shifted  slightly  due  to  an 
intervening  time  interval  or  baseline 
separation 
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A = active  area  of  the  detector 

x and  y = the  horizontal  and  vertical  directions 
The  solution  of  this  integral  for  each  possible  shift  (Ax,  Ay)  gives 
one  point  of  the  correlation  function.  The  particular  values  of 
(Ax,  Ay)  which  give  a maximum  value  for  the  correlation  function 
correspond  to  the  image  misregistrations  (Ref.  8:11.4). 

Recall  that  the  digitized  target  image  lias  been  stored  line  by 
line  in  a memory  unit.  Two  such  memories  are  road  out  with  a delay  in 
one  or  the  other  memory.  The  corresponding  elements  are  multiplied  and 
summed.  The  resulting  number  is  called  the  covariance  function.  The 
maximum  value  of  the  covariance  function  represents  the  highest  degree 
of  correlation  between  the  image  scenes  (see  Figure  8) . The  shift 
which  corresponds  to  the  maximum  value  of  the  covariance  function  is 
the  image  misregistration.  The  maximum  value  is  found  in  one  direction, 
say  the  x direction,  then  the  image  misregistration  is  found  in  the 
other  direction  by  shifting  perpendicular  to  that  point. 

Several  |x">ints  are  computed  in  each  direction.  These  points  can 
bo  used  to  approximate  a correlation  curve  (see  Figure  0).  The  (xiints 
are  run  through  a correlation  algorithm  which  determines  the  exact 
displacement  between  the  two  video  images  (Ref.  8 : 2 . 27)  . 

For  angular  tracking,  the  two  images  compared  are  successive 
images  of  the  target  separated  by  a time  interval.  Difference  in  the 
two  images  is  assumed  to  be  due  to  target  motion  (Ref.  8:2.27) . For 
ranging,  the  images  from  sensor  1 and  2 are  compared  simultaneously. 

The  difference  in  the  two  images  is  assumed  to  be  due  to  baseline 


separation. 


Thus,  range  can  be  computed  by  setting  the  image 


FIGURE  8.  Correlation  Shift  (Ref.  1:228) 


or.  Curve 


. 

misregistration  equal  to  the  parallax  and  using  Equation  (4)  to  compute 
the  range. 

Since  the  target  is  moving,  the  time  interval  between  successive 
images  is  important.  The  smallest  time  interval  jx>ssible  is  the  image 
integration  period.  For  the  Fairchild  CCD  211,  this  is  2.25  ms.  Signal 
read-out  time  is  less  than  integration  time  so  that  there  is  no  delay 
in  picking  up  the  next  integrated  image  signal.  The  additions  and 
multiplications  needed  to  compute  the  covariance  values  plus  the 
correlation  algorithm  also  take  time  to  perform.  In  order  for  the 
digital  processor  to  keep  up  with  the  image  signal  generated  450  times 
every  second,  the  number  of  photosensitive  elements  is  restricted. 

It  was  found  that  a sub-area  of  the  array  of  8 to  10  elements 
squared  centered  on  the  target  was  sufficient  to  perform  adequate 
correlation  (Ref.  8:2.34).  The  maximum  number  of  picture  elements 
allowed  will  be  determined  by  integration  time  requirements  and 
necessary  update  rates  for  the  data.  Since  the  covariance  is  the  sum 
of  the  cross  multiplication  of  the  picture  data,  an  n x m array 
would  take  (n  x m)  multiplications  and  (n  x m)  sums.  If  the  array 
is  increased  by  one  element  in  either  direction,  say  (m  + 1)  , then 

there  are  (n  x m)  + n multiplications  and  sums  to  be  computed  for  each 
relative  shift.  This  means  that  the  processor  time  will  exceed 
maximum  detector  integration  interval  if  the  sub-array  exceeds  some 
threshold  dimensions  (Ref.  8:30). 

The  data  rate  required  to  maintain  track  on  a target  or  update 
data  to  director  and  weapons  release  systems  are  less  than  450  Hz. 

If  these  rates  are  20  Hz  (40  samples/sec) , then  the  data  being  produced 
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at  450  Hz  can  be  averaged  to  increase  accuracy.  A laboratory 
demonstrator  built  by  CAI  has  shown  resolutions  of  1/400  of  a photo- 
sensitive element  at  20  Hz  low  pass  filter  break  frequency  (Ref.  10:1.3- 
1.4)  . 

Summary 

In  conclusion  for  section  IV,  range  accuracy,  AR  , is  related 
to  angular  resolution  of  the  system.  The  angular  resolution  of  the 
system  depends  on  how  the  video  signal  is  produced  and  processed. 

Using  a CCD  focal  plane  array  and  image  correlation,  CAI  Corporation 

. „ , . . photosensitive  element/400 

of  Bourns  Incorporated  claims  a Aa  = i 

This  degree  of  accuracy  is  due  to  the  high  data  rate  of  the  system  and 

assumes  a low  data  rate  requirement  for  tracker  servomechanism , aircraft 

director  system  and  weapons  release  system. 
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V. 


Baseline  Separation  of  Sensors 


Introduction 

By  Equation  (6)  of  Section  II,  AR  , the  range  measurement 
accuracy,  can  be  improved  by  increasing  the  base  between  the  sensors. 

If  the  system  is  to  be  put  on  fighter  aircraft,  the  leading  edge  of  the 
wings  provides  the  widest  possible  base  plus  a relatively  unobstructed 
field  of  view  in  front  of  and  below  the  aircraft.  A possible  exception 
would  be  the  twin  tail  booms  of  the  McDonnell  Douglas  F-15  (Figure  10) 
and  Fairchild-Republic  A-10  (Figure  11) . Note  that  the  downward  view 
of  the  sensors  will  be  blocked  by  the  wings  in  both  of  these  cases.  As 
the  primary  mission  of  the  A-10  is  ground  attack,  this  could  seriously 
limit  field  of  view  (FOV) . Since  the  primary  mission  of  the  F-15  is 
air-to-air,  the  downward  FOV  is  not  as  important.  Table  IV  shows  a 
comparison  of  wing  spans  and  tail  separation. 

Table  IV.  Comparison  of  Aircraft 


' 

Aircraft 

Wingspan 

Vertical  Tail  Separation 

A-10 

57  ft.  - 6 in. 

18  ft.  - 10  in. 

F-15 

42  ft.  - 9.7  in. 

11  ft.  - 9 in. 

F-16/A/B 

31  ft. 

— 

Figure  12  is  a graph  of  maximum  range  as  a function  of  angular 
accuracy  of  the  stereo  system.  The  equation  used  to  generate  this 
graph  is 


(AR/R  ) B 

R = 

max  Aa 

AR/R  = the  range  accuracy  desired  at  maximum  range 
max 

From  Section  III,  this  was  .011 

The  tarqet  must  cover  a window  of  about  seven  by  seven  picture 
elements,  pixels,  so  that  an  adequate  correlation  can  be  performed 
(Ref.  8:2.34).  This  is  required  so  that  the  covariance  function  can 
be  computed  for  a shift  of  +3  units  either  side  of  zero  shift.  If 
the  target  is  a small  fighter,  it  will  display  a cross  section  as 
small  as  one  square  meter  when  viewed  head-on.  The  target  subtends 
6 x 10  *’  radians  at  89NM  in  one  dimension.  In  the  CCD  array,  the 
pixels  are  spaced  every  18  micrometers  vertically  and  30  micrometers 
horizontally.  In  order  that  the  target  cover  seven  of  the  pixels 

with  30  micrometers  spacing,  the  optical  system  must  have  a focal  length 

-4 

of  32.33  meters.  One  element  would  have  a resolution  4.33  x 10  mrad 
in  the  horizontal.  This  is  marked  on  the  graph  in  Figure  12.  By 
insuring  that  the  sensors  are  spaced  by  30  ft,  the  angular  accuracy 
required  of  the  tracking  system  is  larger  than  one  pixel.  In  this  case, 
no  signal  filtering  would  be  required  to  improve  on  the  image  misregis- 
tration found  by  correlation.  Using  a siqnal  processing  to  lower 
Au  to  10  ^ mrad  would  allow  a baseline  separation  of  5 ft  to  range 
accurately  out  to  89NM. 
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Baseline  separations  of  a foot  or  more  will  require  that  the 


ramjinq  sensors  be  housed  in  separate  devices.  In  the  unit  proposed 
by  CA1  , the  ranqe  sensors  optical  axes  were  maintained  parallel  by 
means  of  a sol id  mount.  For  split  sensors,  some  other  method  must  be 
found  to  keep  their  optical  axes  aligned.  The  alignment  must  be  known 
to  within  the  accuracy  required  for  ranging,  Aa  . Unfortunately, 
state  of  the  art  alignment  tor  remotely  located  optical  systems  on  an 
aircraft  is  only  accurate  to  .1  mrad  (Ref  11:8).  Determination  of 
a scheme  to  maintain  optical  alignment  beyond  ,1  mrad  accuracy  is 
beyond  the  scope  of  this  thesis.  Prom  Table  111,  longest  base  available 
is  the  A-10's  ‘30  ft  wingspan.  It  is  apparent  from  Figure  12  that 
An  - l mrad  is  required  for  a base  of  50  feet  and  range  of  89NM. 

Summary 

Maximum  range  capability  using  stereometric  ranging  is  related  to 
tlie  baseline  in  a simple  way  (Kquation  (8)).  Using  a focal  plane 
array  and  correlation  technique  described  in  the  CAI  proposal  requires 
an  optical  system  with  a focal  length  of  32-1/3  meters.  Separating  the 
range  sensors  to  gain  a wider  base  brings  out  the  problem  of  optical 
a 1 iq ration t . 

Once  the  ranging  sensors  have  been  mounted  on  the  wings,  they 
will  be  subject  to  vibrations  which  will  affect  their  performance . 

The  next  section  considers  the  effect  of  structural  vibrations  on 
the  sensors. 
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VI . Effect  of  Flexure  and  Vibration  on  Optical  Sensors 


Introduction 

Aircraft  structures,  especially  the  wings  and  tail,  are  subject 
to  a number  of  loads  and  forcing  functions.  As  a result,  they  flex 
and  vibrate.  Anything  attached  to  these  structures  are  also  subject 
to  the  effects  of  flexure  and  vibration. 

Flexure  and  vibration  are  terms  used  to  describe  two  types  of 
elastic  motion.  Flexure  is  the  elastic  bending  of  a structure  due 
to  changes  in  flight  condition  (Ref.  11:4).  They  are  low  frequency 
bending.  For  our  purposes,  flexures  will  have  a frequency  of  5 Hz  or 
less.  Vibrations  are  higher  frequency  elastic  deformations  caused  by 
such  things  as  wind  gusts,  engine  vibrations,  and  flight  control  motion 
(Ref.  11:4).  For  the  purpose  of  this  thesis,  elastic  bending  at 
frequencies  higher  than  5 Hz  is  considered  to  be  vibration. 

Flexures  and  vibrations  may  be  classified  as  linear  or  angular. 

When  they  are  linear,  they  transport  the  sensor  from  one  position 
to  another  in  a straight  line.  When  they  are  angular,  they  rotate 
the  sensor  about  an  axis . 

In  the  analysis  of  structural  elastic  motion  acting  on  the  sensors, 
a coordinate  system  relative  to  the  sensors  and  target  will  be  used. 

The  plane  containing  the  two  sensors  and  the  target  is  the  sensor-target 
plane.  'Die  direction  normal  to  this  plane  is  the  sensor-target  normal. 
The  direction  from  sensor  to  target  is  called  LOS  for  line  of  sight. 

The  last  direction  is  perpendicular  to  the  other  two  directions  and 
is  called  LOS  perpendicular.  Angular  vibrations  about  these  three 
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Sc-nsor  Coordinates 


directions  will  be  called  ' V2  ' ^3  ^seo  Figure  13) . The  riqht 

hand  rule  applies. 

When  analyzing  these  elastic  motions,  the  effect  of  their  amplitude 
and  rate  on  the  sensor  will  be  considered.  Vibrational  acceleration 
affects  the  sensor  mainly  through  moments  and  torques  on  the  seivo- 
mochanisms.  Analysis  of  the  servomechanisms  is  not  within  the  scope 
of  this  thesis.  The  sensors  are  assumed  to  bo  qimbul  mounted  in 
elevation  and  azimuth  on  the  aircraft. 

The  purpose  of  this  section  is  to  determine  t he  etfect  of  flexure 
and  vibration  on  focal  plane  sensors.  Thi s effect  will  be  expressed  as 
an  equation  whenever  possible. 

Linear  Vibrations  Effect 

Linear  vibrational  effects  are  considered  tirst.  Their  amplitude 
and  rate  in  each  of  the  three  directions  result  in  an  interaction  with 
the  system  which  is  unique  for  each  direction. 

Sensor-Target  Normal . In  the  sensor-target  normal  direction, 
vibrational  amplitude  is  measured  as  changes  in  target  elevation.  Both 
sensors,  when  they  are  tracking,  point  directly  at  the  target.  After 
every  integration  period,  the  sensor  is  corrected  via  the  servomechanisms 
for  any  image  misalignment  found  through  signal  processing.  The 
vibrating  sensor  indicates  that  it  is  the  target  which  is  vibrating  uj 
and  down  about  an  equilibrium  position  (see  Figure  14).  Since  the  sensor 
measures  angles,  a displacement  in  the  sensor-target  normal  is  measured 
as 


(9) 
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where 


$ = angle  measured  about  the  LOS  perpendicular 

P 

A = sensor  displacement  along  the  sensor-target  normal 
N 

R = range  to  the  target  from  the  sensor  of  interest. 

This  motion  has  no  effect  on  ranging.  The  actual  elevation  angle 
is  found  by  averaging  the  elevation  signal  from  the  sensor  over  an 
interval  of  time.  The  averaging  can  be  done  in  the  following  manner. 
Let 

0 = elevation  angle  of  the  target  using  the  baseline  as  an 

axis  (Pig.  14a) 

T = elevation  measured  by  the  sensor 
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if  ij>  (t)  is  a periodic  function  with  a frequency  u>^  > — 
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since  0 (t  ) - constant  over  the  period*!’ 
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where 


T = an  interval  of  time  such  that  — < T < — — ■ 

w V 

v 0 


V,  = rate  at  which  0 changes 

0 p 

P 


iov  = frequency  of  the  vibi'ation 


By  averaging  the  signal,  vibrations  which  go  through  several 

oscillations  in  time  T have  an  average  value  close  to  zero.  For 

2 

the  purpose  of  this  thesis,  vibrations  with  a frequency  of  = — 
are  assumed  to  be  averaged  to  zero.  Vibrations  below  frequency  u>v 
will  be  averaged  together  to  determine  the  mean  amplitude  and  rate 
for  the  vibrating  structure. 

The  vibration  can  be  thought  of  as  a composite  of  vibrations 
with  a continuous  distribution  of  frequencies.  In  theory,  this 
distribution  of  frequencies  goes  from  0 to  . In  reality, 

however,  the  data  collected  on  these  vibrations  only  cover  the  range 
5 to  1000  Hz.  The  lower  limit  is  due  to  the  inability  of  the  sensors 
to  measure  vibrations  accurately  below  that  frequency.  The  upper 
limit  is  observed  empirically  from  the  data. 

The  maximum  value  for  T is  important,  since  it  determines  u 
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The  maximum  value  for  T is  limited  by  the  assumption  that  0 is 

a constant  during  T . In  Table  III,  the  LOS  angular  accuracy, 

, is  specified  to  be  a maximum  of  1 mrad  „ 0^^  , the  LOS 

angle  to  the  target,  is  a combination  of  the  target's  elevation  and 

azimuth  angles.  If  the  azimuth  angle  is  suppressed  for  this  analysis, 

G = 0 . Likewise,  A0  ~ A0  _ and  0=0  . If  0 is 

LOS  p p LOS  p LOS  p 

is  the  rate  at  which  the  target's  elevation  angle  changes,  then 
T = AO  /6  . The  rate  at  which  0 changes  depends  on  several  factors, 

p P P 

such  as  the  type  of  target,  its  velocity,  and  its  range.  Of  the  two 
basic  types  of  missions,  air-to-air  and  air-to-ground,  the  air-to-ground 
mission  usually  involves  the  highest  rates  for  0^ 

For  example,  if  a fighter  makes  a level  bomb  run  on  a ground  target 
(Fig.  14b) 


6 = V.  sin  0 

P f P 


Assume  sin  0 

P 


A 

R 


V 

then  0 = — - 

P R 


where 

= velocity  of  the  fighter 
A = fighter's  altitude  above  the  target 
R = slant  range  to  the  target 
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T 


AO 


AO  R2 


(Ob) 


Figure  14c  is  a graph  of  Equation  (9b)  . The  two  curves  plotted 

in  Figure  14c  represent  two  different  fighter  velocities.  The 

approximate  release  ranges  based  on  bomb  time  of  fall  and  fighter 

altitude  are  marked  by  vertical  lines.  The  value  of  T corresponding 

2 

to  these  ranges  are  marked  with  horizontal  lines.  Since  w = — , 

v T 

w = 71  llz  for  R 7000  ft  and  w = 95  Hz  for  R 5000  ft.  A is 
v v n 

found  by  averaging  the  amplitudes  of  the  vibrations  from  5 Hz  to 

u>  . Therefore,  it  is  necessary  to  locate  the  sensors  at  a position 
v 

where  A is  no  greater  than 


A 

n 

R 

n 


AO 

P 


(10) 


where 

R = minimum  range  capability  of  the  tracker, 
n 

Minimum  range  is  used  to  calculate  A because  it  puts  the  greatest 

n 

restriction  on  A 

n 

Vibrational  rate  lias  another  effect  on  sensor  accuracy.  During 
the  integration  period,  the  target  is  smeared  over  the  detector  array. 

The  amount  of  smear  is  equal  to  the  vibration  rate  times  the  integration 
period  times  the  focal  length,  all  divided  by  the  range.  Image  smear 
affects  sensor  |>ointinq  accuracy  by  distorting  the  video  signal.  The 
relationship  between  image  smear  and  pointing  accuracy  could  be 
determined  by  experimentation.  The  experiment  would  consist  of  increasing 
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FIGURE  14c.  T Vs.  Range  for  Bomb  Run 


the  rate  of  vibration  while  observing  minting  accuracy. 


Given  a 


required  point inq  accuracy,  the  threshold  for  image  smear  would  be  set. 

if  this  threshold  is  represented  by  o , then 

LT 


o 


LT 


H 

n 


t .F 
1 


where 

A = maximum  allowable  linear  vibration  rate  along  the 
sensor-target  normal 

t ^ = integration  period  for  detector 

LOG  Perpend icular . Linear  vibrations  in  this  direction  affect 
each  sensor's  measurement  of  the  tarqet's  azimuth.  Amplitude  of  the 
vibration  is  measured  as  an  angle,  just  as  it  is  for  the  sensor-target 
normal  direction.  Similar  to  Equation  (9) 

A 

e = -£ 

a R (Kef.  Fig.  19) 


where 

o = angle  measured  about  the  sensor-target  normal 
a 

A = amplitude  of  vibration  in  the  LOS  perpendicular  direction 
P 

Similar  to  Equation  (9a),  the  true  azimuth  angle  could  be  found  by 

averaging  the  signal  over  a period  of  time,  T . Since  u>v  = — * 

the  vibrations  above  m will  not  affect  the  measurement  of  0 , 

v a 

the  azimuth  angle  of  tile  target  measured  from  the  fighter's  longitudinal 

is . 


When  measuring  0 , the  maximum  value  of  T is  limited  by  the 

a 

assumption  that  0 is  a constant  during  T . If  it  is  desired  to 

measure  0 to  1 mrad  (Table  III)  and  0 is  known,  then 

a a 


0^  , just  like  0 , depends  on  many  factors  such  as  type  of  target 
and  mission.  Of  the  two  basic  missions,  air-to-air  and  air-to-ground, 
air-to-air  missions  usually  produce  the  greatest  values  for  0 

a 

For  example,  if  a fighter  performs  an  intercept  on  another  aircraft, 
the  0^  can  reach  values  which  significantly  limit  T . An  intercept 
can  be  broken  into  three  phases,  based  on  range.  The  initial  phase 
of  the  intercept  extends  from  initial  range  to  about  10  NM.  This  phase 
of  the  intercept  is  characterized  by  very  low  0 . This  is  because 

cl 

t!ie  fighter  maintains  the  target  at  the  0 whicli  would  cause  them  to 

cl 

eventually  collide.  This  0 does  not  change  at  all  as  the  two  aircraft 
approach  each  other.  The  conversion  phase  of  the  intercept  occurs  in 
the  ten  to  two  nautical  mile  range.  During  this  phase  of  the  intercept, 
the  fighter  maneuvers  from  a head-on  attack  to  a stern  attack.  This 
phase  of  the  intercept  is  characterized  by  relatively  higli  values  for 
0^  . The  final  phase  of  the  intercept  occurs  from  two  nautical  miles 

down  to  minimum  range.  The  fighter  chases  the  target  down  during  this 
phase.  The  angular  l, OS  rate  is  very  low  here  also,  and  his  little 
effect  on  T 

During  the  conversion  phase  of  the  intercept,  the  fighter  usually 


executes  a pure  pursuit  curve.  If  this  is  the  case, 


00 


V 


o = — 

a R 


where 


V 


V = velocity  of  the  target  perpendicular  to  the  LOS 
R = range  to  the  target 
for  a pure  pursuit  curve  is  given  by 


V = V sin  0 

tg  t asp 


where 


V 

tgt 

— 

target's  velocity 

0 

asp 

= 

target's  aspect  angle. 

The  aspect  angle  is  an  angle  determined  by  the  target's  longitudinal 
axis  and  the  LOS  from  the  fighter  to  the  target  (Ref.  Fig.  15a).  The 
aspect  angle  increases  from  0°  to  180°  as  the  fighter  converts 
to  the  target's  stern.  Again,  assuming  the  fighter  executes  a pure 
pursuit  path  as  it  converts  to  the  stern,  0 may  be  approximated  by 


0 

asp 


AO 

a 

Since  T = — 


then 


T 


tgt 


. 1 mrad  R 
sin[~l(R-2)  ] 


2 NM  < R < 10  NM 


A plot  of  this  function  shows  that  T “ .04  seconds  from  3 to  6 
nautical  miles  (see  Fig.  15b).  Therefore,  vibrations  of  frequencies 


(12b) 
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FIGURE  15a.  Aspect  Angle 


Similar  to 


5 to  50  Hz  will  contribute  to  the  mean  value  for 
Equation  (10) , 


A 

-li  = 4>  = AO  '■  (12c) 

K a a 

n 

The  amplitude  of  vibrations  in  the  LOS  perpendicular  direction 
also  affects  the  stereometric  tracker's  range  measurement.  The 
tracker  measures  target  range  by  computing  the  difference  in  target 
azimuth  measured  by  the  two  sensors  (Pig.  4)  . ^ ia  measured  by  each 

d 

sensor  as  [>art  of  the  azimuth  angle  to  the  target.  Figure  16  represents 
two  sensors  mounted  on  a wing  which  is  subject  to  vibrations.  Sensor  2 
is  allowed  to  do  all  the  vibrating  while  sensor  1 is  steady.  In  effect, 
is  the  sum  of  the  angular  deflection  due  to  both  sensors.  ■5'  is 
the  difference  between  0^  and  0 , (see  Fig.  4) . From  Fig.  16 
it  can  be  seen  that 


4*  = (a  - if  ) 


R = R(4>)  = R(a  - <|>  ) 

cl 


R 

n 


ctn 


The  last  equality  uses  the  small  angle  approximation. 


R 


» 


B 

(«-$  ) 
U 


(13) 
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Aircraft 


Since  the  sensors  are  mounted  on  an  aircraft  moving  through  the 
air,  a is  a function  of  time.  It  is  a slowly  varying  function  of 
time  compared  to  <t>  . An  average  or  expected  value  of  R , E{h)  , 

is  given  by 


frequency  at  which  the  sensor  is  vibrating 
rate  at  which  a is  changing 


H 0 * 

fco 

/t  +T' 


* U') 

’ (t')  (1  + rpr— ) df 

a a ( t ) 


a(t'>  dt  + T 


M 


(t')  dt' 


From  the  assumption  on  T , R (t)  and  a(t)  can  be  considered  to  be 

a 

constants  over  the  integration. 


t +T 4 

R itj  r co 


= R (t  ) + 
a 0 


i W f 

T'  »(t0)  J 


4>  (O  dt  ' 
a 


■VT' 


4>  (t')  dt ' = 0 

a 


Since  4>  is  a periodic  function  with  frequency  > T 


(Ref.  12) 
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The  upper  limit  on  the  period  of  time,  T , over  which  e{r} 
is  averaged  is  set  by  the  assumption  that  R^(t)  and  a(t)  are 
constants.  In  actuality,  ci(t)  has  a rate  which  can  be  expressed  in 
terms  of  range , R , as 


_B_ 

o 

R 


V 

c 


where 


• da 

a ~ dt 

V = rate  at  which  the  target  and  fighter  are  closing 
c 

The  accuracy  to  which  a is  known,  Aa  , can  also  be  expressed  in 
terms  of  range: 


Aa 


‘f> 


B 

R 


C 


B 

R 


where 


R = desired  range  accuracy 

(— ) = a constant  for  all  ranges 

R 

The  upper  limit  for  T is  found  from 


T 


A a 
a 


(16b) 


If  the  target  is  another  aircraft  with  a velocity  of  500  kts  and  the 
fighter,  also  at  500  kts,  is  attacking  head-on,  Vc  is  1000  kts. 

If  AR/R  =1.1  NM/R9  NM  = .0125  as  made  a constant  for  all  ranges. 
Equation  (16b)  becomes 
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T 


(16c) 


Figure  15a  is  a graph  of  this  relationship.  For  a ground  target 


the  fighter's  velocity.  Equation  (16c)  becomes  T = (.09 


The  linear  relationship  between  T and  R assumes  the  fighter 


maintains  the  same  aspect  to  the  target  during  the  attack.  If  the 


fighter  is  attacking  another  aircraft,  the  crew  may  want  to  convert 


from  a frontal  attack  to  a stern  attack.  For  the  example  of  a fighter 


attacking  a target  head-on  with  a V of  1000  kts , the  crew  should 


start  a conversion  to  the  target's  stern  around  ten  nautical  miles 


The  crew  will  try  to  maintain  V 


target's  stern  by  2 NM  and  then  continue  to  close  to  minimum  range 


Using  the  relationship  in  Equation  (16b) 


(16d) 


For  the  case  in  which  the  fighter  performs  a frontal  attack  to  a missile 


launch,  the  crew  may  want  the  ability  to  track  accurately  down  to 


minimum  launch  range  of  the  missile.  If  the  minimum  range  is  approxi 


mated  by  10%  of  the  maximum  range,  then  Equation  (16c)  gives 
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FIGURE  16a.  T Vs.  Range  for  Constant  Azimuth 


S i IHV 


lj  Il2  for  a frontal  attack. 


In  cither 


case,  the  amplitudes  of  the  vibrations  from  0 to  to  will  contribute 

v 

the  mean  value  of  A . Thu  position  at  which  the  sensor  is  placed 
will  have  to  be  such  that 


J1  • 

R 


. ,'V<  ", 

A“  1 K It’ 


(17) 


Comparing  the  graphs  ot  figure  lf>b  and  16a  shows  that  T is 
limited  by  llquat ion  ( 16b)  down  to  ranges  of  10  NM.  If  the  crew  converts 
to  the  stern  using  a pursuit  curve,  liquation  (12b)  yields  t he  qreater 
limitation  on  T . Uni  inn  the  final  phase  of  a stern  attack, 
liquation  ( Kid)  is  valid.  These  equations  will  be  used  in  Section  Vll 
to  determine  the  upper  trequeney  limit  of  the  vibrations  which  will 
influence  the  performance  of  the  optical  trackei . however,  it  one 
considers  that  the  sensors  are  held  on  tarqet  by  rate  servomechanisms, 
it  is  joss  ib le  to  use  even  lonqer  inteqration  time  intervals  than 


computed  in  these  eiiuat  ions.  The  rates  0 


0 and  a , are 
P 


relatively  constant.  The  rate  at  which  the  servos  correct  the 

sensor’s  optical  axis  to  the  tarqet’ s IXfS  compensates  for  them.  The 

sensors  respond  to  angular  accelerations,  then,  not  rates.  Since  the 

acceleration  for  the  examples  qiven  are  very  low,  T would  be  qreatly 

increased.  This,  in  turn,  would  bring  <*>v  below  the  frequencies 

figured  previously.  Since  liquations  (9a)  and  (16a)  are  valid  for 

T io  l , vibrations  with  frequencies  from  0 to  w are  used 
v v 


to  compute  A^  and  A 
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rc^apppiS*" 


Vibration  rate  in  tho  LOS  perpendicular  direction  also  causes 
image  smear.  Analysis  in  the  preceding  section  is  valid  here  alos. 
The  equivalent  of  Equation  (11)  is 


% 

A 

P 


o 


LT 


R 

n 


t.F 

1 


(10) 


where 

A = maximum  allowable  linear  vibration  rate  along  the 
LOS  perpendicular 


LOS . The  amplitude  of  linear  motion  in  the  LOS  direction  has  a 
slight  magnification  effect  on  the  target  (see  Fig.  17).  The  area 
which  the  target's  image  covers  in  the  focal  plane  is  expanded  according 
to  the  ratio 


R + Af 
R 


a2  ai  * Aa! 


(19) 


where 


a 

L 


area  of  image  at  time  T^ 
a2  al 

area  of  image  at  time  T 

iunplitude  of  the  vibration  (-)  sign  means  that  A^ 
is  towards  the  target  / (+)  sign  means  that  A 

L 

is  away  from  the  target 


Equation  (18) 


cun  be  written 


R 


6a. 


(20) 


Since  Aa  causes  no  change  in  the  LOS,  vibrations  in  the  direction  of 
LOS  do  not  affect  0 or  range  measurement. 

LlO  b 

The  image  is  smeared  over  an  area  according  to 


a 


1 


AA 


L 


R 


(21) 
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whore 


A A 


L 


t. 

1 


and 

Aj  = rate  of  l.OS  vibration 

t ^ = integration  period 


When  Aa  is  large  enough,  it  affects  correlation  in  the  same  manner 
as  image  smear  did  in  the  other  directions.  Since  the  image  is 
expanding  in  both  directions  at  once,  the  threshold  value  is  expressed 
as  an  area.  The  equation  for  maximum  vibration  rate  takes  the  form 


o _R_ 

a t . 
1 


(22) 


where 

o = threshold  area  referenced  to  a 
a 

a = area  of  the  image  a equilibrium 

A ng ular  Vibra t ion  l: f fects 

Angular  vibrations  are  handled  in  the  same  format  as  the  linear 

vibrations.  Their  amplitude  and  rates  about  each  of  the  three  axes 

are  analyzed  with  respect  to  sensor  performance. 

Angular  Vibration  About  Sensor-Target  Normal.  Angular  vibrations 

about  the  sensor-target  normal  are  v^  . The  amplitude  of  these 

vibrations  directly  affects  the  measurement  of  azimuth  and  range 

(see  Fig.  18).  This  angle  is  in  the  same  plane  as  | of  Equation  (11). 

a 


G4 


The  difference  here  is  that  there  is  no  dependence  on  ranqe  to  the 

target.  Analysis  of  the  motion's  effect  on  the  stereometric  system 

is  exactly  the  same  as  for  *{>  at  frequencies  above  w . Signal 

a v 

averaging  as  in  Equation  (1*1)  does  not  average  out  at  frequencies 

below  to  . Error  in  measuring  u is  equal  to  the  sum  of  the  squares 
of  the  amplitude  of  vibration  for  both  sensors.  Putting  this  in  the 
form  of  an  equation  (see  Fig.  14) 


Act  = /2 


(23) 


where  was  assumed  to  be  the  same  for  each  sensor.  calculated 

from  Equation  (21)  is  maximum  allowable  amplitude  for  angular  vibrations. 

The  rate  of  these  angular  vibrations  multiplied  by  the  integration 
period  and  focal  length  gives  the  length  of  image  smear  in  the  EOS 
perpendicular  direction.  This  can  be  compared  to  the  threshold  value 
for  acceptable  pointing  accuracy.  Solving  the  equation  for  angular 
rate  gives  the  following  equation. 


v 


1 


F t. 

l 


(24) 


where 

= maximum  allowable  vibration  rate. 

Angular  Vibration  About  LOS  Perpendicular.  Angular  vibrations 
about  the  LOS  perpendicular  are  . Amplitude  of  these  vibrations 

has  an  effect  on  the  elevation  measurements  by  the  sensors.  When  these 
vibrations  have  a frequency  greater  than  w 
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, the  vibration  amplitude 


VI 


can  be  averaged  out.  At  and  below,  this  vibrational  amplitude  will 

be  measured  as  pitch  variation.  The  desired  pitch  angle  accuracy, 

AG  , is  set  equal  to  v .In  the  form  of  an  equation,  the  last 
P 2 

sentence  is 


are  superimposed.  Angular  vibrations  about  this  axis  will  make  the 
target's  image  appear  to  rotate  (see  Pig.  19). 


Figure  19.  Image  Rotation 

Analysis  of  the  effect  of  this  kind  of  rotational  smear  on  targets 
of  various  shapes  can  be  done  experimentally.  The  experiment  should 
determine  the  amount  of  image  rotation  allowable  before  correlation 
processing  is  seriously  degraded.  If  this  value  is  designated  as 

°LOS  ' th°n 


' LOS 


(27) 
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r 


where 

threshold  angle  at  which  correlation  processing  is 
degraded 

maximum  allowable  angular  rate  of  vibration  about  LOS 
axis 

This  concludes  the  analysis  of  vibrations.  Flexures  are  analyzed 
in  the  next  section. 

Flexures 

Flexures  are  the  same  movements  as  vibration,  only  at  a lower 
frequency.  For  this  system,  vibrations  were  considered  low  frequency 
at  iuv  or  less.  Flexures  fall  into  this  category;  therefore,  they 
have  already  been  analyzed  in  the  vibration  section. 

Summary 

This  section  analyzed  the  effects  of  vibrations  and  flexures  on 

sensor  performance.  The  goal  of  this  analysis  was  to  find  parametric 

expressions  for  the  effect  of  these  vibrations  and  flexures.  A 

coordinate  system  was  established.  Linear  and  angular  vibrations 

were  distinguished  from  each  other. 

Vibrations  were  analyzed  first.  Linear  vibrations  along  the 

sensor-target  normal  affected  elevation  tracking  performance.  Averaging 

the  sensor  data  will  eliminate  vibrations  at  frequencies  greater  than 

u>  . Limiting  amplitude  for  vibrations  below  w and  flexures  is 
v v 

given  by  Equation  (10)  . 


A = R [A0  ] 
n n p 


(10) 


69 


The  limiting  rate  is  calculated  from  Equation  (11) 


o R 
LT  n 


t.  F 
1 


Linear  vibrations  along  the  LOS  perpendicular  affect  azimuth  and 

range  measurements.  Averaging  sensor  data  will  eliminate  vibrations 

above  . For  vibrations  below  and  flexures, 

v v 


R AO  , for  azimuth 
n a 


<f> 

v=r  R 


Vibration  rate  is  set  by  Equation  (18) 


, for  range 


P 


o R 
LT  jn 

t . F 

l 


Linear  vibrations  along  the  LOS  have  a slight  magnifying  effect. 
The  rate  of  vibration  is  limited  by  Equation  (21) . 


o R 

a n 

a t . 

l 


(11) 


(12a) 


(17) 


(18) 


(22) 


Angular  vibrations  about  the  sensor- target  normal  affect  ranging. 

Averaging  sensor  data  will  eliminate  vibrations  above  m . Vibrations 

below  u are  flexures  are  limited  in  amplitude  by 
v 


Aa_ 

J2 


(23) 
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Maximum  angular  rate  is  given  by 


v 


1 


F t. 
1 


(24) 


Angular  vibrations  about  the  LOS  perpendicular  affect  elevation. 
Amplitude  for  vibrations  below  and  for  flexures  is  given  by 

Equation  (24).  Maximum  angle  rate  is  given  by  Equation  (25). 
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AO 
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(25) 
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F t . 
1 


(26) 


Angular  vibrations  about  the  LOS  direction  do  not  affect  range 
or  elevation  tracking  directly.  They  can  cause  image  smear.  Equation 
(26)  gives  the  maximum  angular  rate  for  vibrations  or  flexures. 


v 


3 


o 

LOS 

t . 

l 


(27) 


Since  the  main  distinction  between  flexures  and  vibrations  is 
frequency,  the  equations  derived  for  vibrations  are  valid  for 
flexure . 

In  the  next  section,  these  equations  will  be  used  to  aid  in  the 
placement  of  sensors  on  a tactical  aircraft. 
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VII.  Analysis  Applied  to  F-15  Data 

In Product  ion 

This  section  takes  the  analytical  expressions  of  Section  V]  and 
uses  them  to  determine  optimum  sensor  placement.  The  equations  in 
Section  VI  determine  maximum  allowable  vibration  amplitudes  and  rates. 
Expressions  in  this  section  determine  vibration  amplitude  and  rate  of 
the  sensor  as  a function  of  position  when  possible.  Therefore,  it  is 
[lossible  to  find  a position  of  the  sensor  such  that  amplitudes  and 
rates  are  less  than  the  maximum  allowable. 

The  leading  edges  of  the  wings  are  the  primary  candidates  for 
location  of  the  sensors.  The  twin  tail  structure  of  the  F-15  is  also 
considered.  An  exact  relationship  between  vibrations  and  position 
along  the  leading  edge  of  the  wing  would  require  analysis  by  moans 
of  finite  elements.  Analysis  of  this  type  had  not  been  accomplished 
on  a fighter  aircraft  previously,  and  such  an  analysis  is  not  within 
the  scope  of  this  thesis. 

Some  flight  tent  data  is  available  for  F-15  wing  and  tail 
structures.  The  flight  tests  were  conducted  to  determine  the  maximum 
vibrational  environments  at  certain  positions  on  the  wing.  The  data 
collected  on  these  flight  tests  is  in  the  form  of  power  spectral  density, 
r:q>  (Figs.  20,  21,  22).  In  order  to  find  the  mean  value  for  amplitude 
in. I r.ite  el  f he  vibrations,  it  is  necessary  to  use  the  relationships 
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I.  I I V.  I III  Appendix  A: 


r “> 


(27a) 


V (t) 


I / *. 

Lo  w1-  J 


[G] 


X (t) 


uo 

*2  / 
-0 


pa(w)  1* 


4 

u) 


[G] 


(27b) 


where 


V (t) 
X (t) 
P (w) 

cl 

G 


= _ average  velocity  of  the  vibrating  sensor 
= average  displacement  of  the  vibrating  sensor 
= power  spectral  density 

= earth's  gravitational  acceleration  constant  at 
sea  level 


The  factor  of  G is  required  because  the  data  was  originally  recorded 
in  G's  squared. 

All  the  data  available  is  from  linear  accelerometers.  Angular 
vibration  data  is  not  available  because  adequate  angular  accelerometers 
are  not  available  (Ref.  17) . 

Data  collected  for  the  F-15  wing  and  tail  represent  severe 
vibrational  environments.  Data  representing  normal  flight  conditions 
is  better  suited  to  this  analysis.  However,  such  data  is  not  available 
at  this  time.  The  rates  and  amplitudes  computed  from  the  available  data 
should  be  considered  as  maximum  values  and  not  the  normal  values.  Two 
types  of  maneuvers  were  performed  to  generate  the  vibrations.  They 
were  a windup  turn  and  a symmetric  pullup.  Of  these  maneuvers,  the 
windup  turn  induces  the  greatest  vibrations.  When  possible  then,  the 
data  for  the  symmetrical  pullup  maneuver  will  be  used  in  preference  to 
the  windup  turn  data. 
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FIGURE  20  Power  Spectral  Density  for  LU10 
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The  limit  of  integration  for  Equations  (27 a)  and  (27b)  represent 
the  entire  spectrum  of  frequencies.  The  data,  however,  is  limited  to 
a minimum  of  5 Hz  because  the  accelerometers  used  to  measure  vibrational 
accelerations  had  a minimum  frequency  range  of  5 Hz.  However,  it  is 
likely  that  there  are  no  vibrational  modes  of  the  wing  or  tail  below 
5 Hz. 

If  the  wing  is  modeled  as  a tapered  cantilevered  beam,  the 
frequency  of  the  first  mode  of  vibration  is  given  by  (Ref.  16:456-470) 


f 


b 

1 


(~) 

3p 


l. 

2 


where 

b = width  of  the  beam  at  the  support 
1 = length  of  the  beam 

E = modulus  of  elasticity  of  aluminum 

p = density  of  the  beam  density  of  aluminum 

Density  of  the  beam  is  assumed  to  be  approximately  one-tenth  of  the 
density  of  solid  aluminum.  Width  of  the  wing  at  the  wing  root  is 
approximately  8 inches  when  viewing  the  wing  head-on.  Width  of  the 
wing  at  the  root  is  225  inches  when  viewing  the  wing  in  planform  as 
in  Figure  23. 

Width  of  the  tail  at  the  fuselage  is  6.5  inches  when  viewed  head-on. 
The  wing  is  256  inches  long  from  fuselage  to  wing  tip.  The  tail  is 
124  inches  high  from  the  fuselage.  Using  these  figures,  the  frequencies 


of  the  first  mode  of  vibration  are 
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FIGURE  23.  F-15  Wing  Top  Planform  View 


wing  z 


150  Hz 


f . 

wing  x 

f . , 
tail  y 


18  Hz 


where 

f . = frequency  of  wing  vibrations  in  the  z direction 

wing  z 

f . = frequency  of  wing  vibrations  in  the  x direction 

wing  x 

f = frequency  of  tail  vibrations  in  the  y direction 

tail  y 

The  frequency  computed  for  the  first  mode  of  wing  vibration  in  the 

x direction  is  obviously  too  high.  This  is  due  to  the  inaccuracy  of 

the  model  used  for  wing  vibrations  in  this  direction,  f . also 

tail  y 

seems  high.  However,  f . and  f . , are  certainly  greater 

wing  x tail  y 

than  f . . The  cantilevered  beam  model  for  the  wing  vibrating 

wing  z 

in  the  z direction  is  probably  a good  approximation  because  b is 

so  much  smaller  than  1 . The  estimations  on  p and  E will  not 

change  f . more  than  an  order  of  magnitude.  For  the  purpose  of 

wing  z 

this  analysis,  then,  the  minimum  frequency  for  wing  or  tail  vibrations 
in  any  of  the  directions  considered  is  5 Hz. 

The  upper  limit  on  frequency  of  vibration  is  . In  Section  VI , 

wv  is  shown  to  have  a range  of  values,  depending  on  T . T , in 
turn,  depended  on  the  type  of  target  and  attack.  From  Figures  14c, 

15b,  and  16a,  jT  is  an  interval  of  time  between  4 seconds  and  .021 
seconds.  In  cc*nputing  X and  V from  the  data,  = 95  Hz  gives  the 

largest  values,  while  = .25  Hz  yields  the  least.  Since  the  frequency 


of  the  first  vibration  mode  of  any  of  the  structures  is  5 Hz,  an 


u>  below  that  frequency  yields  X and  V of  0 The  two  extremes 

for  make  analysis  of  the  data  somewhat  uncertain.  A more  reasonable 

v 

cut-off  for  is  10  Ilz.  This  frequency  represents  a period,  T , 

of  .2  seconds  which  is  possible  for  most  parts  of  the  missions  discussed 
in  Section  VI. 

Table  V is  a list  of  the  values  for  X and  V calculated  at 
positions  LU  10  , LU  11  , LU  33  and  the  tail.  See  Figures  23  and 
24  for  these  locations  on  the  aircraft.  The  accelerometers  at  positions 
LU  10  and  LU  11  measured  acceleration  in  the  z direction  (Fig.  23). 

The  accelerometer  at  position  LU  33  measured  accelerations  in  the  x 
direction  and  the  tail  accelerometer  measured  accelerations  in  the  y 


direction.  X and  V were  computed  using  symmetric  pullup  data  at 
positions  LU  10  and  LU  11  , while  windup  turn  data  were  used  at  the 
other  two  locations. 

If  the  values  for  X and  V at  w = .25  Hz  are  designated 

v 

as  a , and  at  - 10  Hz  as  m , and  at  95  Hz  as  b , then  the 


of  the  sensors. 


In  analyzing  the  effects  of  vibrations  on  the  sensors,  it  is 

convenient  to  clearly  describe  0 and  0 . These  angles  are 

pa 

depicted  in  Figure  25.  Azimuth,  0 , is  measured  in  the  x-y  plant 

cl 

using  a line  in  the  x-y  plane  perpendicular  to  the  azimuth  line  and 
through  the  origin. 

The  following  subsections  use  the  flight  test  data  and  equations 
from  Section  VI  to  determine  suitable  placement  of  the  electro-optic 
sensors.  Analysis  will  start  with  the  twin  tails. 

Tail  Vibrations 

Sensors  located  at  the  position  of  accelerometer  in  Figure  24 

undergo  linear  vibrations  of  15  inches  in  the  y direction.  Assuming 

both  tail  structures  behave  the  same,  each  sensor  will  undergo 

displacements  of  (15  inches)  cos  0 = A 

A p 

Equation  (17)  gives  the  relationship 


A 

— = AO 
R a 

n 


where  AO  = AO  is  given  as  1 mrad  in  Table  III.  If  A =15  inches, 

a LOS  p 

then  R = 1250  ft.  This,  then,  is  the  minimum  range  of  the  stereometric 
n 

tracker  mounted  on  the  tail. 

The  direction  of  the  amplitudes  of  these  vibrations  is  along  y , 


which  is  also  the  direction  of  the  base  line.  Therefore,  A = AB 

P 


If  the  errors  in  range  due  to  a and  F are  ignored,  then 


FIGURE  25.  Aircraft  Coordinate  System 


From  Table  III,  R = 1.1  NM  when  R = 89  NM  . Solving  Equation  (28) 

for  A gives  A =6  inches  if  B = 40  ft. 

P ' P 

If  a is  equal  to  the  horizontal  dimension  of  one  pixel,  the 
LT  1 

maximum  allowable  rate  of  vibration  is 

• o R 

A = — — £ .41  ft/sec  (18) 

p t F 

Image  smear,  then,  will  cause  the  greatest  problem.  In  view  of  the 
27  ft/sec  rate  from  the  data  for  the  tail  pods,  they  seem  to  be  an 
unsuitable  place  for  the  sensors. 

Wing  Vibrations 

A search  of  F-15,  F-16  and  A-10  flight  test  data  produced  limited 
data  for  linear  vibrations  of  the  F-15  wing.  RMS  amplitudes  and  rates 
of  vibration  are  given  in  Table  V.  Linear  vibration  data  is  recorded 
in  the  x and  z direction. 

In  order  to  have  a relationship  between  vibrational  amplitude  and 
position  along  the  leading  edge  of  the  wing,  it  was  necessary  to  model 
the  wing  as  a cantilevered  beam  of  variable  cross  section  (Ref  16:465-470) . 
The  details  of  the  analysis  can  be  found  in  Appendix  B.  By  means  of 
this  analysis,  the  first  mode  of  vibration  is  approximated  by 

Vq  = <29> 


where 

d = distance  from  wing  tip  to  position  of  interest 
1 = length  of  the  wing 
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= amplitude  coefficient 
y^  = perpendicular  amplitude  of  the  beam 


The  results  of  this  analysis  must  be  used  with  caution.  The 

model  falls  short  of  the  wing  in  a number  of  ways.  First,  the  wing 

vibrates  in  three  dimensions,  not  just  two  independent  planes.  The 

wing  was  modeled  as  a symmetric  wedge.  In  reality,  it  is  a complex 

form  with  little  symmetry  in  the  x - y plane.  As  can  be  seen  from 

the  PSD  charts,  there  are  many  vibrational  frequencies.  RMS  vibration 

amplitude  is  computed  from  the  data  and  taken  as  the  amplitude  of  only 

the  first  mode  of  the  bending  beam.  This  will  give  too  large  an 

amplitude  for  the  first  mode.  The  model  is  justified  because  it  is 

the  best  available  to  determine  a relationship  between  d and  y^ 

given  the  lack  of  available  data  for  the  wings.  Witn  these  things  in 

mind,  linear  and  angular  vibrations  will  be  analyzed. 

The  wing  will  be  analyzed  at  position  LU  10  first.  Linear 

vibrations  below  w in  the  z direction  cause  errors  in  measuring 
v 

the  elevation  angle.  Fran  Table  III,  the  accuracy  desired  is  .1  mrad. 
A^  = 3.2  inches  from  Table  V.  Using  Equation  (10) , 


R 

n 


[A6  ] 
P 


2.66  ft 


A can  be  increased  6 inches  if  an  R = 500  ft  is  desired. 

N n 

Using  Equation  (29)  where  it  is  understood  that  y1  is  in  the 
z direction,  a^  can  be  found  by  letting  y^  = An  and  d = the 
position  of  LU  10  . z is  the  amplitude  of  the  vibration. 
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4.6  inches 


(29a) 


31  (1  - d/l)2 


If  Equation  (29)  is  solved  for  d and  if  y^  is  set  equal  to  the 
desired  amplitude  of  vibration,  d will  represent  the  minimum  distance 
from  the  winq  tip.  Thus,  Equation  (29)  becomes 


d = 


(30) 


If  y = A^  , then  d = 29  inches,  which  is  147  inches  from  the  fuselage. 
Since  the  fuselage  is  160  inches  wide,  B = 38  ft 

If  is  equal  to  the  vertical  dimension  of  one  pixel,  the 

maximum  vibration  rate  will  be 

• °LT 

A = — — = .67  ft/sec  (11) 

N t . 

1 

where 

A = Z cos  0 
N p 

Z = mean  of  the  vibrations  in  the  z direction 

Again,  this  is  a severe  restriction  since  Z was  3.5  ft/sec. 

Unfortunately,  no  equations  like  (29)  or  (30)  exist  for  vibrational 
rates.  If  they  did,  it  would  be  possible  to  determine  d such  that 
the  vibration  rate  in  the  z direction  is  acceptable. 

Vibrations  in  the  z direction  will  give  rise  to  A as  the 

Li 

sensor  follows  a target  in  elevation. 


(31) 


A (sin  0 ) (cos  0 ) 
a p a 


When  the  sensor  is  looking  directly  up  or  down,  0 = 90°  ,0=0, 

p a 

and  A . Vibration  in  this  direction  has  no  effect  on  tracker 

t.  a 

per  fortnance . 

The  other  direction  of  linear  wing  vibration  is  in  the  x direction. 
Vibration  in  the  x direction  is  measured  by  the  accelerometer  at 
1,0  33  . Using  the  RMS  value  of  x , amplitude  of  the  vibration  in 
tlie  x direction,  as  y ^ , a = 3.3  inches  by  Equation  (29). 

The  amplitudes  and  rates  of  wing  vibrations  are  related  to  sensor 
vibrations  by 


A - X sin  0 
N a p 


A,,  = X sin  0 
N a p 


A - X s in  0 
P a a 


A = X sin  0 
P a a 


A X cos  0 cos  0 
L a p a 


A = X cos  d cos  0 
L a p a 


where 


rate  of  vibration  in  the  x direction 


The  requirements  on  the  amplitudes  of  sensor  vibration  remain  the 

samej  it  is  only  the  direction  of  vibration  which  has  changed.  Thus, 

• • 

the  requirements  A = 6 inches  , A = .67  ft/sec  , and  A = .41  ft/sec 
n n P 

remain  unchanged.  For  a baseline  of  40  ft,  A =5.9  inches 

P 

From  Table  V,  X .13  inches  . The  location  of  the  accelerometer 

a 


which  produced  this  data  was  115  inches  from  the  aircraft  centerline. 


If  a baseline  of  <10  ft  is  desired,  the  mean  .;implitude  of  the  vibration 

at  the  new  position  can  be  computed  using  Equation  (29a).  First, 

solving  for  a,  at  LU  33  (d  = 141  inclies,  y,  = X ) 
l la 

X 

a 

a 3.M  inches 

(1  - d/l) ^ 

if  the  two  eleetro-opt ic  sensors  are  placed  40  ft  apart,  the  mean 
value  of  X will  be  increased  to 


X = a (1  - d/1) ‘ = 2.6  inches 

a 1 


Structural  vibrations  of  this  amplitude  are  within  the  limit  set  by 


A and  A 
P n 

There  are  no  angular  vibration  data  for  the  wing.  Equation  (29) 
can  be  used  to  transform  the  linear  data  to  angular  data  in  the 
following  manner: 
s ince 

yx  = f(d) 

where  y^  and  d are  orthogonal 


0 = 


d(d) 


d (d) 


f(d) 


-2  a 


i ( i - f i 


H9 


(35) 


r 


where 

0 = angle  of  the  tangent  at  the  point  d measured  in  radians 

For  example,  the  angular  deflection  at  d *»  l.U  10  , in  the  y - z 

plane,  is  0 = .13  rad  = 7.25"  . 

Amplitudes  and  rates  of  wing  angular  vibrations  are  related  to 
sensor  angular  vibrat ions  by 

— — • • • 

v,  =0  sin  0 cos  0 , v.  * 0 sin  0 cos  0 (36) 

1  z p a 1 z p a 


— • • • 

v_  = 0 sin  0 cos  0 , * 0 sin  0 cos  0 (37) 

2  z a p 2 z a p 


v,  = 0 cos  0 cos  0 , v_  = 0 cos  1)  cos  6 (38) 

3  z p a 3 z p a 


From  Section  VI , angular  vibrations  with  frequencies  below 


% 


have  a maximum  amplitude  of  v = — - (Eg.  (23)) , v = 5.6  x 10  J rad  , 

1 /a  1 

and  for  v>  = AO  (Eq.  (25))  , v„  = . 1 mrad  . The  angul  - rates  of 

2 •’  o 2 „ 

those  vibrations  are  * — — (Eq.  (24)),  = 2.47  x 10  rad/sec  , 

ci  i o 

• LT  . -a  . LT 

v..  = — — (Eq.  (26)),  v = 4.12  x 10  rad/sec,  v = — — (Eq.  (27)), 


F t. 

i 

= 9.07  rad/sec 


t. 

i 


if  the  electro-optic  sensors  have  a baseline  of  40  ft, 


0^  = 4.7  x 10  rad  . By  Equation  (36) 


0 sin  0 
z p 


v when  the 
-2 


target  is  directly  in  front  of  the  fighter.  Since  0^  «=  4.7  x 10 

-5  — 

and  v,  =*  5.6  x 10  rad  , an  elevation  angle  oi  0 * arc  sin  v,/0  = 

1 p 1 z 

1.2  x 10  rad  is  the  maximum  allowable  before  0 begins  to  influence 

z 


the  range  accuracy  of  the  tracker.  If  0 were  as  high  as  30°  , the 

P 

position  of  the  sensors  would  have  to  be  adjusted  so  that  0 ■ 2v^ 


90 


From  equation  (35)  and  0 = 2 

■ 


1 


<t:>(2v  ) 

d = — - l = -176  inches  (39) 

2ai 

The  minus  sign  has  no  significance  here.  In  order  to  meet  the  require- 
ment of  , the  sensor  must  be  placed  next  to  the  fuselage. 

Since  the  cantilevered  beam  model  does  not  transform  angular  rates, 
there  are  no  comparisons  to  make . 

For  wing  vibrations  in  the  x - v plane,  0 = 7.6  x 10  ^ = .43° 

x 

using  Equation  (35)  and  d = l.U  33  . 0 is  related  to  the  sensor 

x 

amplitudes  by 

v,  = 0 cos  0 (40) 

1 x p 


V-  = 0 sin  0 (41) 

3 x p 

As  in  the  case  for  0 , 0 can  be  equal  to  v,  or  \>  given  the 

z x 13 

proper  azimuth  or  elevation  angle.  Even  though  0^  < 0^  , it  still 

is  not  small  enough  to  meet  the  restriction  of  . Using  Equation 

(39) , d a 176  inches  , which  is  next  to  the  fuselage.  However,  0^ 
is  more  restrictive,  so  that  the  value  of  d it  predicts  must  be 
honored . 

Summary 

Tills  section  used  the  equations  in  Section  VI  and  data  from  flight 
tests  to  determine  optimum  positions  for  the  optical  trackers.  The  data 
available  was  in  the  form  of  linear  PSD.  A simplified  program  calculated 


RMS  amplitudes  and  rates  from  the  data.  It  was  necessary  to  model  the 
wing  as  a cantilevered  beam  in  order  to  determine  sensor  placement. 

Analysis  of  the  tail  data  indicated  that  the  amplitude  of  the  linear 
vibrations  is  too  large  by  a factor  of  15:0  . The  mean  rate  of 

vibration  at  the  tail  was  also  larger  than  that  allowed  for  an  image 
smeat  of  one  pixel.  The  tail  is  not  a suitable  location  for  the  electro- 
optic  sensors  because  of  the  vibration. 

Wing  vibrations  wore  analyzed  with  the  aid  of  Equation  (30). 
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(30) 


The  mean  amplitude  of  linear  vibration  in  either  the  x - y or  y - z 
plane  is  :.mu  1 1 enough  that  it  does  not  affect  the  tracker  performance. 
The  moan  angular  rate  (1.0  10)  is  an  order  of  magnitude  higher  than 

A 

n 

Equation  (3M)  was  used  to  find  an  acceptable  position  in  terms  of 
ur.gulai  vibration  amplitude. 
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Vj  was  the  most  restrictive  anqular  amplitude.  It  caused  the  sensors 
to  ho  placed  next  to  the  aircraft  fuselage. 

It  must  be  remembered  tluxt  the  flight  conditions  under  which  this 
data  wore  taken  are  worst-case  vibrational  environments.  Normal  flight 
conditions  would  cause  much  less  vibration.  Also,  the  model  used  to 


let  ermine  placement  of  the  sensors  on  the  wing  is  severely  limited  in 


VII  l . Conclusions  and  Recommendat ions 


Summary 

The  objective  of  this  thesis  is  a generalized  parumetr ic  trade-off 
analysis  of  a passive  stereometric  tracker.  The  stereometric  tracker 
uses  two  focal  plane  array  sensors  to  track  the  target  LOS.  Azimuth 
and  elevation  are  measured  directly  from  the  sensor's  LOS.  Range  is 
computed  from  Equation  (4)  , R = B/a  , where  a = 0 ^ - 0^ 

Analysis  of  tracker  performance  centered  on  the  system's  angular 
resolution,  baseline  separation  and  the  effects  of  aircraft  vibrations 
and  flexures  on  range  finding  accuracy.  Performance  criteria  were 
specified  in  Section  III  and  were  transformed  into  sensor  criteria  in 
Section  Vl.  The  effects  of  aircraft  vibrations  and  flexures  on  tracking 
performance  were  analyzed  in  Section  VII  and  put  into  analytic  expressions. 
Flight  test  data  was  used  to  find  optimum  placement  of  optical  sensors 
on  an  F-15. 

In  Section  IV,  correlation  between  the  sensors  was  found  to  be  the 
most  accurate  means  of  determining  a . Once  the  sensors  are  split, 
however,  a means  of  maintaining  optical  alignment  to  1 mrad  is 
necessary.  The  means  of  measuring  alignment  to  1 mrad  between  two 
electro-optic  sensors  were  not  investigated. 

Assuming  that  some  means  is  found  either  to  measure  or  maintain 
optical  alignment  between  sensors,  range  tracking  can  be  improved  by 
increasing  baseline  separation  (Ref.  Equation  (6)).  Aircraft  structures 
best  suited  are  the  wings  and  the  twin  tail  design  of  some  fighters. 

Those  structures  are  subject  to  severe  vibrational  and  flexure  environments. 
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Very  little  data  are  available  on  these  structures  to  characterize  their 
response  to  flight  conditions.  Information  which  is  available  is  in  the 
form  of  linear  PSD.  Angular  vibration  information  is  critical  in 
evaluating  the  performance  of  an  optical  system.  Some  evaluation  can 
be  made  using  techniques  described  in  References  17,  18,  19  and  20. 
However,  as  stressed  in  those  references,  these  are  not  substitutes 
for  direct  measurement  of  angular  vibration  data. 

Dy  modeling  the  aircraft  wing  as  a tapered,  cantilevered  beam, 
a trade-off  between  baseline  separation  and  angular  vibration  can  be 
made.  Optimum  position  of  the  sensor  using  these  two  factors  requires 
that  the  sensors  be  placed  next  to  the  fuselage.  This  is  where  A a 
is  as  small  as  required  for  the  desired  range  accuracy.  However,  this 
position  is  unsuitable  due  to  the  interference  of  the  fuselage  with 
the  tracker's  field  of  regard. 

Conclusions 

1.  The  means  of  maintaining  or  measuring  alignment  between  remotely 
placed  optical  sensors  is  a major  problem  to  be  overcome . 

2.  A more  accurate  relationship  between  position  and  vibration 
at  that  position  is  required. 

3.  Angular  vibration  data  is  required  to  adequately  analyze  the 
system.  Angular  vibrations  were  much  more  detrimental  to  the  tracker's 
performance  than  the  linear  vibrations. 

4.  Vibrations  and  flexures  below  u drive  the  design  of  the 

v 

system  to  a single  unit  with  the  optical  sensors  rigidly  mounted. 

5.  The  passive  stereometric  has  the  potential  to  give  an  aircraft 
the  element  of  surprise  when  attacking  a target.  The  element  of  surprise 
is  particularly  useful  for  tactical  aircraft. 
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Kocommend.it  ions 

The  following  items  should  be  considered  before  evaluting  remotely 
{laced  sensors  for  stereometric  tracking. 

1.  A means  of  measuring  or  maintaining  optical  alignment  to 
1 mrad  should  be  investigated. 

2.  A finite  element  analysis  should  be  made  on  the  wing  and  tail 
structures  to  compute  angular  as  well  as  linear  vibration  data. 

3.  A means  of  collecting  angular  vibrations  of  aircraft  structures 
in  flight  should  be  found  to  validate  recommendation  2. 

4.  The  investigation  into  extending  the  range  of  a passive 
stereometric  tracker  should  be  continued  because  of  its  tactical  advantages 
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Appendix  A 

I n teq ration  of  Power  Spectral  Density  (Ref . 12) 


Flight  test  data  is  collected  in  the  form  of  acceleration 
history  on  tape.  The  acceleration  data  is  transformed  into  power 
spectral  density  by 
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To  compute  the  velocity  from  acceleration  data 
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Appendix  b 

F irst  Mode  of  Vi  I iration  For  a Cantilevered  boom 
(Ref.  lb : 4 56-470) 


The  following  derivation  of  the  first  mode  shape  of  a cantilevered 
beam  can  be  found  in  the  reference  cited. 

The  differential  equation  for  a vibrating  beam  in  the  shape  of 
a wedge  has  the  form  (see  Figure  2b) 
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= first  moment  of  inertia 
area  of  the  cross  section 


length  of  the  beam 
depth  at  the  fixed  end 
flecural  rigidity  of  the  beam 
mass/unit  area 


The  boundary  conditions  are 
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Fitjure  26.  Beam  (Kef.  16:466) 


The  following  series  satisfies  these  boundary  conditions  using  the 
Ritz  Method  (Ref  16:461) . 


y = 


As  a first  approximation 


Kurt  F.  Schroeder  was  born  on  10  January  1947  in  Sparta, 
Wisconsin.  After  graduating  from  high  school  in  Seoul,  Korea  in  1900, 
he  returned  to  the  United  States  to  attend  Lamar  University  where  he 
received  a Bachelor  of  Science  degree  in  Physics  in  May  1970.  Upon 
completion  of  OTS  in  San  Antonio,  Texas,  he  was  commissioned  in  the 
USAF.  He  completed  flight  training  and  received  his  wings  in  1971. 

He  served  as  weapon  systems  operator  in  the  F-4  in  the  13  TPS,  Udorn, 
Thailand,  and  in  the  490th  TPS,  Hahn,  Germany  until  entering  the 
School  of  Engineering,  Air  Force  Institute  of  Technology  in  June  1977. 


Permanent  Address : 1444  Avenue  E 

Nederland,  Texas  77627 


unclassified 

SFCUNITV  CLASSIFICATION  OF  THIS  PACE  7>*Ti»,.  limit  fnlrlmj) 


REPORT  DOCUMENTATION  PAGE 


I REPORT  number  2 GOVT  ACCESSION  NO 

AFIT/GEP/PH/78D-11 


4 TlTLf  nod  Subtitle) 

PARAMETRIC  ANALYSIS  OF  STEREOMETRIC  TRACKER 
FOR  USE  IN  TACTICAL  AIRCRAFT 


I 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOG  NUMBER 


AU  THOR(  »; 

Kurt  F.  Schroeder 
Capt  USAF 


9 PERFORMING  ORGANISATION  NAMl  ANO  ADORE  SS 

Air  Force  Institute  of  Technology  (AFIT/EN) 
Wr ight-Patterson  AFB  OH  *15433 


ii  controlling  office  name  ano  aooress 

Air  Force  Avionics  Laboratory  (AFAL/RWT) 
Wright-Patterson  AFB  OH  45433 


14  MONITORING  AGENCY  NAME  a AOOHfSS<  it  ditterent  from  Controlling  Othc*) 


5 TYPE  OP  REPORT  A PtRlOO  COVEREO 

MS  Thesis 


6 PERFORMING  ORG.  REPORT  NUMBER 


8 CONTRACT  OR  GRANT  NUM0CR(») 


10  PROGRAM  ELEMENT.  PROJECT.  T ASK 
ART  A A WORK  UNIT  NUMBERS 


99912004 


12  REPORT  OATE 

October  1978 


IV  NUMBER  OF  PAGES 

104 


15  SECURITY  CLASS.  ( of  (Mi  report) 


Unclassified 

15*~d1eCL  A SSI  F*  I CAT  1 0 N OOWNG  R A Dm  G 
schedule 


IS  DISTRIBUTION  STATEMENT  (ol  thin  Report) 

Approved  for  public  release;  distribution  unlimited. 


»8  SUF*F*{_  EMEN  T AR>  NOTES 

19  KE  Y WOROS  fConllniia  on  reverae  aide  H necehmmrv  end  Mertfffy  by  block  number) 

Optical  tracking 

Image  degradation 

Correlation 

Passive  Range  Finder 

\i 

;0  ATTRACT  (Conrlnu#  »w*  reverse  aid*  it  nece«*«rv  end  identity  hr  bfock  number) 

The  stereometric  ranqe  finding  technique  is  used  as  the  basis  of  an  air- 
borne passive  tracker.  An  analysis  of  the  method  reveals  that  angular  accuracy 
is  the  most  critical  element  in  this  technique.  Range  accuracy  can  be  improved 
by  increasing  baseline  separation.  However,  the  vibration  environment  of  the 
aircraft  becomes  worse  as  the  baseline  is  increased.  The  results  of  this 
investigation  indicate  that  these  vibrations  counter  any  gain  in  resolution. 

The  investigation  was  hindered  by  a lack  of  vibrational  data,  especially 


I A 71  '473  EOITION  OF  lltov  »»  IS  OBSOLETE 


UNCLASSIFIED 

"'ll  T > • I l I * • 


